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ABSTRACT 



Context. This work is conducted as part of the "SPECPDR" program dedicated to the study of very small particles and chemistry in photo- 
dissociation regions with the Spitzer Space Telescope (SST). 

Aims. We study the evolution of the Aromatic Infrared Bands (AIBs) emitters across the illuminated edge of the Horsehead nebula and 
especially their survival and properties in the HII region. 

Methods. We present spectral mapping observations taken with the Infrared Spectrograph (IRS) at wavelengths 5.2-38 jim. The spectra have a 
resolving power of /1/A/l = 64- 128 and show the main aromatic bands, H2 rotational lines, ionised gas lines and continuum. The maps have an 
angular resolution of 3. 6-10.6" and allow us to study the nebula, from the HII diffuse region in front of the nebula to the inner dense region. 
Results. A strong AIB at 11.3 fim is detected in the HII region, relative to the other AIBs at 6.2, 7.7 and 8.6 fim, and up to an angular separation 
of ~ 20 "(or 0.04 pc) from the ionisation front. The intensity of this band appears to be correlated with the intensity of the [Nell] at 12.8 fim and 
of Hot, which shows that the emitters of the 11.3 fim band are located in the ionised gas. The survival of AIB emitters in the HII region could 
be due to the moderate intensity of the radiation field (Go ~ 100) and the lack of photons with energy above ~25 eV. The enhancement of the 
intensity of the 1 1.3 yum band in the HII region, relative to the other AIBs can be explained by the presence of neutral PAHs. 
Conclusions. Our observations highlight a transition region between ionised and neutral PAHs observed with ideal conditions in our Galaxy. A 
scenario where PAHs can survive in HII regions and be significantly neutral could explain the detection of a prominent 11.3 fim band in other 
Spitzer observations. 

Key words. ISMiindividual objects: IC434, Horsehead - ISM:dust, extinction - ISM: HII region - Infrared: ISM - ISM: fines and bands 



1. Introduction 

Polycyclic Aromatic 
posed twenty ye ars 



Hydroc arbons (PAHs) were 



ago by iLeger & Puget 



pro- 
(11984 and 



AUamandola et al.l (119851) to explain the infrared emission 
bands observed at 3.3, 6.2, 7.7, 8.6 and 1 1.3 fim. These emitters 
are an ubiquitous component of interstellar dust which has 
been observed with the Infrared Space Observatory (IS O) in 



a wide range of inte r stella r conditions (e. g. Boulanger et al 



1998al: lUchidaetall ExXll) . These aromatic infrared bands 
(AIBs) have already been observed in spectra attributed 



to HII regions (e.g . iPeeters et al.L l2002t iHonv et all 12001 



Vermeil et al.L 120021) but never with clear proof that their emit- 



ters are within the ionised gas rather than within an associated 
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Telescope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under a contract with NASA. 



photodissociation region on the same line of sight. Moreover, 
several studies report the destruction of these emitters in the 
HII regions of M17 and the Orion Bar. In these regions, the 
strong radiation field is thought to b e the main cause of this 
destruction (e. g. iKassis et al. [ 120061 and reference therein). 
In this paper, we use "PAHs" as a generic term in order to 
designate the emitters of the AIBs although the exact nature of 
these emitters is still a matter of debate. 

Dust grains can play an important role in the ener- 
getic balance of HII regions through photoelectric heating 
( Weingartner & Draine , 200ll and references therein). Since 
these processes are dominated by small grains, the presence 
of PAHs in HII regions has a strong impact on the physics of 
these objects. 

In front of the western illuminated edge of the molec- 
ular cloud L1630, the visible plates are dominated by ex- 
tended red emission due to the Hg line emission emerging 
from the HII region IC434 (e. g. iLouisel 119821) . In the visi- 
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Fig. 1. Ha map, in arbitrary units, obtained with t he 0.9 m Kitt 



Peak National Observatory (KPNO) telescope dPound et al 



2003h . The contours show the areas observed with both IRS- 
SL (small area) and IRS-LL (large area). 




Fig. 2. Average SL spectra obtained within the Inner and HII 
regions (see the contours of these two regions in Fig.|5]) before 
(thick line) and after (thin line) correction of the residual emis- 
sion (see §|2]l. 



ble, the Horsehead nebula, also known as B33 (lBarnardLll919l) . 
emerges from the edge of L1630 as a dark cloud in the near 
side of IC434. The Horsehead nebula is a familiar object 
in astronomy and has been observed many times at visible 



IR and submm wavelengths (IZhou et al 



2003 



2005 



Pound et all 20031: Tevssier et al . , ., 

Petv et al.Ll2005HHilv-Blant et alil2005h . IC434 and the 



1993; Abergel et al. 
12004: Habart et al. 



Horsehead nebula are ex cited by the cr Or i onis s tar which is 
an 09.5V binary system ( Warren & Hesse4 19771) with an ef - 
fective temperature of ~ 34 600 K (.Schaerer & de Koteilll997b . 



L1630 is located at a distance of ~400 pcQ. Assuming that cr 
Orionis and the Horsehead are in the same plane perpendicu- 
lar to the line of sight, the distance between them is ~3.5pc 
(~0.5°) which gives Go ~ 100 (energy density of the radiation 
field between 6 and 13.6 eV in unit of Habing field, iHabind, 
1968h for the radiation field which illuminates the Horsehead 
nebula. 

In this paper, we study the AIBs from IC434 in front of 
the Horsehead nebula observed with the Infrared Spectrograph 
(IRS; H ouck et al., 2004) on board the Spitzer Space Telescope 
( Werner et al ■i l2004h . The paper is organised as follows : in §|2] 
we present our IRS data and the data processing. In §|3] we 
show that IRS data confirm the description of the structure of 
the object from previous studies. In the following section (§|4|i, 
we extract the typical spectrum of the HII region and study the 
location of the AIB emitters of this spectrum. We compare the 
HII region spectrum with the spectrum obtained in the inner 
region in §|5] and propose a scenario to explain the observed 
spectral variations in §|6] The survival of PAHs in the HII re- 
gion is discussed in §|7] We conclude in §|8] 

2. Observations & data reduction 

The Horsehead nebula has bee n observed w i th IRS as a part 
of our "SPECPDR" program djobhn et all l2005b on 2004 



October 2, and using the Short-High (SH), Long-High (LH) 
, Short-Low (SL) and Long-Low (LL) modules of the in- 
strument. In this paper, we only present SL (5.2-14.5;um, 
slit size: 57" X 3.6", R=64-128) and LL (14-38//m, slit size: 
168" X 10.6", R=64-128) observations. We used the "spectral 
mapping mode". An observation is made of Nstep= 23 (SL) or 
25 (LL) steps of half the slit width in the direction perpendicu- 
lar to the slit long axis. For the SL module, three observations 
were taken successively at three different positions on the sky 
in order to perform a complete mapping of the illuminated edge 
of the nebula. The resulting observed areas are shown in Fig.[T] 
overplotted on the Ha map. The integration times were 14 and 
60 s per pointing for the second (5.2-8.7 fim) and the first (7.4- 
14.5;um) orders of SL, respectively, and 14 s per pointing for 
both orders of LL. 

We have developed a pipeline which builds spectral cubes 
(two spatial dimensions and one spectral dimension) in a ho- 
mogeneous way from the data (version SI 3) delivered by the 
Spitzer Science Center (SSC). We start from the BCD level. 
One integration corresponds to one BCD image. For each BCD 
image, we extract for all wavelengths an image of the slit which 
is projected on the sky. For each observation, made of Nstep in- 
tegrations, we build a spectral cube with x Ny spatial pixels 
and Nw spectral pixels. We keep the same wavelength sampling 
as in the BCD images and the spatial grid has a pixel size of 
2.5" for LL and 0.9" for SL (which corresponds to half the 
pixel size on the BCD images). Whenever we study the full 
spectral range, we also reproject the SL data on the LL grid. 

We identify and correct the bad pixels not flagged out 
in the SSC pipeline by median filtering on the combined 



' from the study of the d istances to B stars in the Orion association 
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Fig. 3. Maps of the Horsehead nebula at 11-11.5 (SL, left panel) and 20-2 1 jim (LL, right panel). The SL and LL maps have 
pixel size 0.9"and 2.5", respectively. Positions of the spectra of Fig.|4]are shown on the LL map (right panel). 



Nstep BCD images. The data are flux-calibrated in Jansky us- 
ing the pipeline S 1 3 conversion factors and the tuning coeffi- 
cients given in the "fluxcon" table. Finally, we derive extended 
emission flux intensities by using the Slit Loss Correction 
Functio n due to the point-spread function overfilling the IRS 
slit (Smi th et all 12007). In the following, the LL data at 35- 
38 jum are not considered due to the strong decrease of sensi- 
tivity. 

Only the HII region in front the Horsehead nebula presents 
nearly flat infrared emission and could have been used to de- 
fine an spectrum. Since the goal of this paper is precisely to 
study the emission emerging from the HII region, we did not 
subtract any off spectrum. The lack of such correction explains 
the discontinuity systematically found for all pixels between 
the first order and the second order parts of the SL spectra, 
together with a systematic decrease of the continuum with de- 
creasing wavelength, down to negative values for wavelengths 
below ~6 jum (Fig.©. We find that for each pixel the ampli- 
tude of these effects does not depend on the detected emission, 
but appears strongly correlated with the non-zero emission de- 
tected in the interorder regions of the BCD images which in 
principle does not receive any incident photon. This non-zero 
emission does not depend on the emission in the order regions 
of the BCD images, but presents some correlation with the 
emission in the peakup region of the array. It could be a residual 
after the "droop" or the stray-light corrections (see the IRS data 
handboolJl). We estimate the amplitude of this residual in the 
order regions by extrapolating for each row the emission de- 
tected in the interorder regions, and subtract this residual from 
the BCD image. This subtraction is performed on the BCD im- 
age multiplied by the flat field image (taken in the calibration 
files delivered by the SSC), since we work in the hypothesis 
that the effects we want to correct are additive. Then, we divide 
the corrected BCD image by the flat field image. Finally, we 
build a corrected spectral cube using the algorithm described 



above. The correction only affects the shape of the continuum 
emission and does not change the amplitude of the spectral 
bands and lines. The discontinuity between the two orders of 
the SL spectra is strongly reduced (Fig.©. However, we have 
to keep in mind that at that time the origin of the corrected 
effects is not known, therefore we must be cautious in the in- 
terpretation of the continuum emission. 



3. The Horsehead as seen by IRS 

Figure [3] gives an example of the spectral maps obtained from 
SL and LL observations. For all pixels within the SL field, we 
have a full spectrum from 5 to 35 yum (Fig.HJ. 

The zodiacal emission is computed using the SSC back- 
grou nd estima toj^ which is based on the COBE/DIRBE model 
(Kels all etall [l998). The dashed lines in Fig.g] show the 
zodiacal contribution to our observations which is not sim- 
ply the zodiacal emission at the date and sky coordinates 
of our observations. Indeed, the "dark" level subtracted from 
the data in the SSC pipeline is obtained without shutter by 
pointing a fixed area of the sky with faint infrared emission 
(RA = 268 ?96, DEC = 65?43) as explained in the IRS Data 
Handbook. This "dark" level will contain some zodiacal emis- 
sion (~ 14MJy sr ' at ~ 18/im). Therefore, the zodiacal contri- 
bution to our observations is the difference between the zodi- 
acal emission at the time and the position of our observations 
and the zodiacal emission at the "dark" position. This zodiacal 
contribution does not vary across the observed area of the sky 
and is accurate to ~I4%^. 

The spatial structure detected with IRS is comparable to the 
broad-band observations taken with ISOCAM (Abergel et al.[ 
20031), but we now have the spectral information from 5 to 
35 ;um and better spatial resolution. Thanks to the edge-on ge- 
ometry of this PDR, it is possible to perform spectral analysis 
of the emission in the HII region, the edge of the PDR and the 



^ see _http://ssc.spitzer.caltech.edu/irs/dh/| 



^ see http://ssc.spitzer.caltech.edu/documents/background/| 
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Fig. 4. Typical spectra obtained by combining SL and LL data for individual pixels (1) at the infrared peak position, (2) in the 
inner region of the Horsehead nebula behind the emission peak and (3) in front of the dense cloud in the HII region, as shown 
in the right panel of Fig. [3] Alternating thicknesses correspond to SL2, SLl, LL2 and LLl orders. The dashed lines show the 
contribution of the zodiacal emission (see §[3]l. 



inner region inside the PDR, separately. Three illustrative spec- 
tra for individual pixels are presented in Fig.lD 



(1) The first spectrum is taken at the infrared peak position and 
is typical for a PDR. It shows the main H2 rotational lines 
(0-0 S(4) to S(0) at 6.9, 9.7, 12.3, 17.0 and 28.2;um), the 
AIBs and continuum. 

(2) The second spectrum is taken in the inner region behind the 
peak (to the east of the peak) and shows AIBs and H2 emis- 
sion lines with lower intensities since the emitting matter is 
located more deeply in the dense cloud. 

(3) The third spectrum is taken in front of the illuminated sur- 
face in the HII region (to the west of the peak) and is 
dominated by fine structure lines of ionised species as ex- 
pected for a HII region. It shows the following lines: [Aril] 
at 6.98 jum, [Nell] at 12.8yum, [SIII] at 18.7 and 33.4jum, 
[Sill] at 34.8 //m, but not the more excited lines: [Nelll] 
at 15.5 jum, [SIV] at 10.5 yum and [Arlll] at 9.0iJ.m. It also 
contains the 11.3 fim AIB. 

Spectra (1) and (2) also present ionised lines which are 
likely emitted by the ionised medium surrounding the dense 
cloud. The continuum emission at wavelengths lower than 
~ 20//m appears to be dominated by the zodiacal emission for 
the spectra (2) and (3). 

In the following, we remove the zodiacal contribution from 
all spectra and focus our study on the SL spectra of the HII 
region which contains the 11.3 /urn AIB. The study of the spec- 
tral properties around the peak position will be the subject of a 
forthcoming paper. 



4. HII region spectrum 

Fig.|5] shows emission maps of the 11.3Ajm AIB (fitte d with 
a lorentzian profile following iBoulanger et al. , Il998bh . H2 0- 
0S(2), [Nell] at 12.8yum (both fitted with a gaussian pro- 
file) and Ha obtaine d at the Kitt Peak National Observatory 
(KPNO) telescope by lPound et all (.2003) (Fig.fB- Both [Nell] 
and Ha are emitted by ionised gas. However, they cannot be 
used alone to define the HII region and exclude PDR emis- 
sion which could be located on the same line of sight behind 
or in front of the ionised gas (with respect to us). We see in 
Fig.|5] that the Ha and [Nell] emissions do not peak at the 
same location due to projection effects and the difference of 
extinction efficiency since they emit at different wavelengths. 
We need a PDR tracer in order to exclude PDR emission. We 
use the H2 0-0S(2) line at 12.3 /vm sin ce it is a good PDR 



tracer for this dense illuminated ridge jHabart et al.L l2005h 



which is well detected in our data. Thus, we define the HII 
region as the area where Ih» > 450 in arbitrary units and I(H2 0- 



S(2)) < 6 10 W m sr (detection limit due to the noise in 
individual spectra). We use Ha rather than [Nell] to trace the 
ionised gas since it is not affected by ionisation fraction effects. 
The contours of the HII region we have defined are shown in 
Fig.E] The Ha emission in this HII region traces the ionised 
material in front of the Horsehead. In fact, the extinction due 
to this material on a line of sight is Ay ~ 0.01-0.06, consid- 
ering ne ~nH ~ 100 - 350 cm"^ (see appendix [A| . a depth of 
~0. 1 pc for the Horsehead ne bula jHabart et all 20051) and 
Nh/Av ~ 1.87 10^1 cm-2mag-i (iBohlin et allll978l) . 

The average spectrum computed within this area (Fig.|6l 
lower spectrum) shows the [Aril] and [Nell] lines at 6.98 fim 
and 12.8 jum, respectively. A band is also detected at 11.3jum 
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Fig. 5. Emission maps of the 11.3jum AIB, H2 0-0S(2), [Nell] at 12.8 yum (as seen by IRS) and Ha JPound et all l2003h . The 
contours show the area defined as the HII region (west of the peak) and the inner region (east of the peak) in §|4]and §|5] 




Fig. 6. Mean spectra of the HII region (lower) and the inner region (upper) as defined in §|4] and §|5] and shown in Fig.|5] The 
[Nell] line at 12.8 fj.m is truncated for clarity. In the HII region spectrum, the broad feature at ~ lOfim is an artefact (see §|4]and 
Fig-I21>- The zodiacal emission has been removed from these spectra (see §[3]l. 



with a surprisingly high intensity compared to the intensities of 
other AIBs at 6.2, 7.7 and 8.6 /im. Other bands seem necessary 
to account for the emission plateau between 11.3 and 13/im. 
A broad emission feature is visible around 10 //m, however its 
amplitude is variable (Fig.|7]i and depends strongly on the off- 
set position from the center of the slit (i.e. on the position on 
the detector). We therefore conclude that this feature is mainly 
due to an artefact. On the other hand, the AIB at 11.3/im is 
real since it appears everywhere in the HII region with a strong 
amplitude, as shown in Fig.|2] 



A method to check whether the 11.3/im emitters are lo- 
cated in the ionised gas is to look for spatial correlation be- 
tween the 11.3/vm AIB and the ionised gas lines emission in 
the HII region. The figure|8]presents the 11.3 yum AIB vs [Nell] 
and vs Ha correlation diagrams for pixels of the HII region. 
The computed [Nell] line intensity is not affected by a possible 
contribution of the 12.7yum AIB, since our spectral decompo- 
sition takes into account this AIB. We see that the 1 1 .3 yum AIB 
intensity is correlated with both the [Nell] and Ha ones which 
shows that the observed correlation is not dominated by sys- 
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Fig. 7. Spectra from three areas (small black rectangles) shown on the 11-11.5 fim map (top, middle and lower panels correspond 
to the north, the central and the south area, respectively). The zodiacal emission has been subtracted (see §[3]l. For clarity, the 
[Nell] line is truncated. The broad feature at ~ 10/vm in the spectrum of the top panel is an artefact (see §|4|l. 



tematic effects in the IRS data. The correlation coefficients are 
0.56 and 0.63 for 1 1.3//m AIB vs [Nell] and 1 l.Bjum AIB vs 
Ha, respectively. Note that we have only considered the points 
within the 3cr limit with respect to the linear fitting (dashed 
lines in Fig.[8l). The relatively low values of the correlation co- 
efficients are due to systematic effects on the IRS computed 
maps (striping. Fig.© and to the fact that the relationship be- 
tween the ionised gas lines intensity and the AIBs intensity is 
not systematically linear 

At least about half the 1 1.3 yum emitters on the line of sight 
should be located in the ionised gas in front of the Horsehead 
nebula where the Ha and [Nell] emissions rise, since the two 
correlations range from ~ 2.5 to ~ 4.5 10"^ W m"^ sr"' for the 
1 1.3 jum A IB intensity. A typi cal Cirrus spectrum has been es- 
timated by lFlagey et aP (l2006h from ISOCAM data (Fig.|9l up- 
per spectrum). It presents a shape comparable to the HII spec- 
trum, except for the 6.2 yum AIB which indicates that the con- 
tribution of Cirrus emission does not dominate the detected 
HII spectrum. In any case, the contribution of Cirrus should 
be comparable for the HII region and inner region spectra (as 
defined in the next section), and cannot affect our conclusions 
based on spectral variations between these two spectra. 

5. Comparison with tlie inner region spectrum 

As for the HII region (§|4|i, we define an inner region where 
Ihq- < 135 in arbitrary units and Iv.ii-ii.s^m < 15 MJy sr"' in or- 
der to avoid the bright infrared filament (Fig. ID were the [Nell] 
emission peaks due to projection effets. The average spectrum 
computed within this area (Fig.|6] upper spectrum) shows the 
0-0S(3)-9.7jum and 0-0S(2)-12.3jum H2 rotational lines, the 
AIBs at 6.2, 7.7, 8.6 and 1 1.3 yum and the [Nell] Hne at 12.8 yum. 
The 0-0 S(5)-6.9 jum H2 and [Aril] (6.98 yum) lines are blended. 

The presence of [Nell] and [Aril] lines indicates that this 
spectrum contains ionised gas emission and could contain AIBs 



emitted in the background and foreground ionised medium. 
However, such a contribution from HII region AIBs will not 
change the conclusions of our analysis of the spectral varia- 
tions that we study in the following. 

The relative intensity of the AIBs presents striking differ- 
ences between the HII region and inner region spectra (Fig.|6]l. 
While the intensity of the 11.3 yum band is comparable, the 6.2, 
7.7 and 8.6/im bands present a spectacular decrease in the HII 
region (more than a factor 2-3). The 6.2 and 7.7//m bands are 
attributed to C-C stretching mode, the 8.6/im band to in-plane 
C-H bending modes, and the 11.3 fj.m band to C-H out-of-plane 
bending mode. In the following section, we study the different 
processes which could explain such spectral variation of the 
6-9 yum/ 1 1.3 yum ratio. 



6. Interpretation of the spectral properties 

6.1. Hydrogenation state 

Hydrogenation state of PAHs can be traced by r elative intensi- 



ties o f the 1 1.3, 12, 12.7 and 13.6yumbands (e. g. lSchutte et al 



1993h which are C-H out-of-plane bending modes with one, 
two, three and four H atoms on the same aromatic rin g, re- 
spectively (e. g. lLeger & Pugetlll984lHonv et afl.l200lb . Both 
inner region and HII region spectra (Fig.|6l) present an 11.3 fim 
band and an emission plateau between 11.3 and 13yum above 
the continuum, which can be attributed to hydrogenated PAHs 
(high hydrogenation coverage of PAHs has already b e en re- 



ported in HII r egions and PDRs, e. g. IVermeii et al.L 12002 



Hony et al.Ll2001i) . Moreover the intensity of the 1 1.3 yum band 



compared to the plateau does not present any strong variation 
between the two spectra, which indicates that the hydrogena- 
tion states are comparable. We conclude that hydrogenation ef- 
fects are likely not the main process which could explain the 
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Fig. 8. Relationships between the [Nell], Ha and 11. 3 //m for 
pixels of the HII region (defined in §|4] and shown in Fig.|5]l. 
Contours are histograms of the point densities. The red straight 
lines are hnear fits to the data. The red dashed lines are the 3cr 
limit with respect to the linear fitting. 

difl'erence in the 6-9 yum / 1 1 .3 //m ratio between the two spec- 
tra. 

6.2. Size distribution 



Using the model of IVerstraete et al.l (l2001h . we find that the 
emission ratio of 6-9 fim/ 11. 3 fim is reduced by a factor of 3 
only if the size distribution contains exclusively PAHs bigger 
than lO"' C atoms (~ 30 A) in the HII region while it is a clas 



<■ 



Z) 
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Fig. 9. Cirrus spectrum observed with ISOCAM dFlagev et al 
.2006) and HII region spectrum as defined in §|4] convolved 
with a gaussian yielding a resolving power of 40, comparable 
to ISOCAM-CVF. The 9-10. 5jum range is suppressed since it 
contain an artefact (§|4]i. 



ratio. Neutral PAHs emit significantly less at 6-9 fim than 
at 11-13 fim with respect to charged ones (both anions and 
cations). The inner region spectrum which comes from neutral 
gas presents a high value of the 6-9 fim / 11.3jum ratio, which 
can be explained by charged PAHs (anions or cations). On the 
contrary, the low value of the 6-9 fim / 1 1 .3 fim ratio in the HII 
region spectrum can be explained by the presence of neutral 
PAHs. The spectra extracted from ISOCAM observations 
(4-16jum) of NGC7023 present comparabl e relative intensit y 

l2005h . 



variations attributed to charge effects ( Rapacioh et al 



Our spectra from the HII region and the inner region (Fig.© 
could therefore be attributed to PAH" and PAH^, respectively. 

The charge state of PAHs is mainly determined by the bal- 
ance b etween photoionisation and re c ombination rates o f elec 



trons dWeingartner & Drainel 120011: iBakes & Tielensl 119941) 



which is generally described by the ratio of the UV inten- 
sity to the electronic density, Go/ng. The presence of posi- 
tively charged PAHs in the inner region can be explained by 
(1) the presence of UV photons which efficiently ionise the 
PAHs and (2) a lack of free electrons for the recombination 
{njnu ~ [C]/[H] ~ lO"'*, C^ being the main provider of elec- 
trons in the PDR). 

For the HII region, we use version 05.07 of Cloudy 



sical size distribution (mean size ~ 6 A, e.g. 
19941) in the inner region. Then, we conclude that a change in 
the size distribution due to destruction of smallest species can 
not explained the 6-9 fim/ 11.3 fim ratio variation. 



6.3. Charge state 



dFerland et al.L 119981) in order to der ive a quantitat i ve es 



Bakes & Tielensl timate of the charge s tate o f PAHs dvan Hoof et all 12004 



Weingartner & Drainel 1200 Ih in a fully (ng ~ nn) ionised 



Theoretical 

in 



(e.g. 



2001a b; Bauschliche^ 
Szczepanski et Vala. il993h works 



Langhofi , 
' 2002h and 



19961: 



Bakes et al 



experimental (e. g. 
show that the charge 
state of PAHs has a strong impact on the 6-9 fim/ 11.3 fim 



medium. We perform a simple model with an incident radi- 
ation field defined fo r an 9.5V star of the Costar catalogue 
dSchaerer & de Koten . Il99'7h and located a t 3.5 pc. The gas is 
taken to be at T = 7500 K dFeriaiidi l2003h and with ne ~ 100- 
350 cm"^ (see appendix lAll. For PAH s with radius from 4.5 t o 
10.5 A and distributed as n(a)oca"^-^ dBakes & Tielensl [l994l) . 
we obtain a mean charge of 0.55 - 0.75 electron per PAH, corre- 
sponding to a fraction of neutral PAHs in the HII region of 25 - 
45%. Moreover, the 6-9 fim/ 11.3 fim ratio of the HII region 
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spectrum is in agreem e nt with those predicted by the emission 
model of iBakes et alj (1200 lab for such a charge distribution. 
We conclude that the HII region spectrum can be explained by 
a mixture of neutral and anionic PAHs. 



7. PAHs survival in the HII region 

The 11.3jum band is observed in the HII region up to a dis- 
tance of ~20"(or 0.04 pc) from the ionisation front (see for 
instance Fig.|5]). This distance can be translated to a lower 
limit of the survival time of the emitters equal to ~5 10"' 
years when considering that the gas just ionised at the ioni- 
sation front expands freely in the HII region at the sound speed 

Cs= ,^^/^~10kms-' for T = 7500 K,ju = 0.7 and 7 = 5/3. 

We have seen in §[3]that in the HII region, all species with 
an ionisation potential (IP) lower or equal than that of [SIII] 
(IP=23.34eV) are detected while species with IP higher or 
equal than that of [Arlll] (IP=27.63 eV) are not detected. Thus, 
the incident radiation field in the HII region contains UV pho- 
tons with energy up to ~ 25 eV. Thus, some aromatic emitters 
can survive to a radiation field with Go ~ 100 and photons up 
to ~ 25 eV. For comparison, in more highly excited HII regions 
such as the Orion Bar, the radiation field is more intense (Go 
larger than lO'*) and also harder since the exciting star is an 
06 (0' OriC, e.g. iKassis etall l2006t lAllers et al.Ll2005b with 
Teff ~ 44000 K. The re, PAHs are destr oyed on time scales 

l2006l) . 



lower than 1000 years ( iKassis et al 



It must be emphasized that the presence of AIB emitters in 
the ionised gas could be related to the continuous injection of 
"fresh" matter due to photoevaporation of the Horsehead. 



contain a significant amount of neutral PAHs. On the con- 
trary, PAHs from the inner region must be positively charged. 
Variations of the size distribution of PAHs could also affect the 



6-9 //m/ 11.3 fim ratio (e. g. IVerstraete et al.L 1200 lb but with a 



lower amplitude than charge state variations. 

In galaxies, the presence of neutral PAHs has been sug- 
gested by Kaneda et al. (2005) to explain the prominent emis- 
sion feature at 11.3//m compared to the 6.2, 7.7 and 8.6/zm 
features observed with Spitzer in several elliptical galaxies. The 
present IRS observations of the Horsehead nebula thus provide 
a textbook example in our Galaxy of the transition region be- 
tween ionised and neutral PAHs and allow to derive a physical 
scenario in order to interpret extra-galactic spectra. 

Appendix A: Determination of the electronic 
density 

The electronic density in the HII region is estimated using 
the intensity ratio of the [SIII] lines at 18.7 and 33.4yum. 
These two lines are fitted with gaussian profiles in or- 
der to compute their integrated intensities for all pixels in 
the HII region (see §|4] and Fig.|5]l. The mean intensities 
are I([SIII], 19 yum) = 2.69 + 0. 1 1 10"^ W m^^ sr"' and I([SIII], 
33yum) = 3.91+0.4210"'^Wm-2sr-'. Error bars include the 
calibration uncertainties (absolute stellar model (lcr = 3%), 
the repeatability (~3%) and the uncertainty of the calibra- 
tion at 33.4jum (-10%)). We have I([SIII], 19jum)/I([SIII], 
33 M m) =0.69 ±0.08. Co nseq uently, using atomic cons tants 
from Galavis et aP (Il995h and lMendoza & Z eint^en (19821) and 



assuming a gas temperature of 7500 K (IFerland i2003i) . we de- 
rive an electronic density Ug ~ 1 00-350 cm~^. 



8. Conclusion 



References 



Our main observational result is the detection of a strong 
1 1.3 //m emission band in the HII region facing the Horsehead 
nebula. The spectral imaging capabilities of IRS allow us to 
show that the integrated intensity of the 1 1.3 fim AIB is corre- 
lated with those of the [Nell] line and Ha in the HII region. 
Moreover the spectral variations of the AIB spectrum is clearly 
spatially correlated with the variations of physical conditions 
between the PDR and the HII region. Thus, to our knowledge, 
this is the first time that we detect the presence of AIB emit- 
ters in ionised gas. Consequently, the 11.3 jjm emitters are not 
efliciently destroyed by the incident UV photons which have 
an energy below ~ 25 eV as suggested by the detected ionised 
species. The survival of PAHs in the HII region could be due 
to the moderate intensity of the radiation field (Go ~ 100) and 
the lack of photons above ~ 25 eV, compared to more highly 
excited HII regions (in terms of intensity and hardness), where 
PAHs can be destro yed on time scales lower than 1000 years 



(IKassis et al.L 120061) . It could also be related to the continuous 



photoevaporation of the Horsehead nebula which bring "fresh" 
matter into the ionised gas. 

The enhancement of the intensity of the 11.3//m band in 
the HII region, relative to the other AIBs, can be explained 
by the presence of neutral PAHs. Our modelling of the charge 
state of PAHs with Cloudy confirms that the HII region should 
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